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Damped mechanical systems with various forms of damping are quantized using the path integral 
formalism. In particular, we obtain the path integral kernel for the linearly damped harmonic 
oscillator and a particle in a uniform gravitational field with linearly or quadratically damped 
motion. In each case, we study the evolution of Gaussian wave packets and discuss the characteristic 
features that help us distinguish between different types of damping. For quadratic damping, we 
show that the action and equation of motion of such a system has a connection with the zero 
dimensional version of a currently popular scalar field theory. Furthermore we demonstrate that the 

, equation of motion (for quadratic damping) can be identified as a geodesic equation in a fictitious 

■ two-dimensional space. 
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I. INTRODUCTION 



The presence of damping in a mechanical system is a natural occurrence. For example, consider a particle falling 
through a fluid under gravity. The form of the damping force depends on the value of the Reynolds number Re = ph/rj, 
where p is the fluid density, rj the viscosity coefficient, £ the characteristic length scale, and v the speed of the particle 
in the fluid. For Re < 1 we may assume linear damping and for 10'^ < Re < 2 x 10^ it is more natural to assume 
• that the damping is quadratic in the velocityii The assumption of a nonlinear damping term makes the equation of 
^ ' motion nonlinear and more difficult to handle in general. Several classical mechanical systems with linear as well as 
0^ ' nonlinear damping are exactly solvable.-^ 

CO ' In contrast, the quantum mechanics of damped mechanical systems is not as easy to understand. Usually we write 
, down Schrodinger's equation for a given potential and obtain the energy eigenvalues and eigenfunctions either exactly 
^ ' or by approximation methods. This procedure does not work for damped systems because of either the explicit time 
, dependence or the complicated form of the Lagrangian and hence the Hamiltonian. 

' The primary goals of this article are to show that Lagrangians can be constructed for simple damped systems, 
\ to use these Lagrangians to construct the path integral kernels for damped systems, and to study the wavepacket 
evolution using these kernels. Our results supplement the existing literature on exact path integrals for mechanical 
systems. 

, Our examples include the damped simple harmonic oscillator and the freely falling particle in a uniform gravitational 
■ field in thepresence of linear/ or quadratic damping. Earlier work on the quantization of damped systems can be found 
^ in Refs. [HjH using a variety of techniques such as variational methods, the Fokker-Planck equation, and canonical 
quantization. Path integral techniques have been used by several authors (see for example, Refs. '^.7'»!9l fl0lfllllT^ll3h . 
^ \ A com preh ensive and up-to-date analysis on various aspects of path integrals (with associated references) is available 
• --j '■ in Ref. lil. 

Kn] ■ In Sec. II we outline the path integral formalism, which we shall use extensively. Then in Sec. IIII Al we consider 
^ ; the damped harmonic oscillator and discuss the construction of the kernel for the under-damped case in detail. As 
■ ■ ■ ' a second example, we consider in Sec. IIIIBI a particle falling under gravity in the presence of a linear damping force. 

In Sec. IIVI we focus on quadratic damping in an analogous way. In all these systems we study wavepacket evolution 
and show how the dispersion of the packet provides us with a way of distinguishing between the magnitude as well 
as various forms of the damping force. As an aside (and a motivation for the reader who wishes to find a taste of 
advanced physics from an elementary standpoint) we connect the quadratic damping scenario with a recently studied 
field theory. In the same spirit, we also illustrate how quadratically damped motion can be viewed as a geodesic 
motion in a fictitious two-dimensional space. In Sec. |V] we conclude with a summary of our results. 



II. PATH INTEGRAL FORMALISM 



Before we begin our discussion of the path integral treatment of damped mechanical systems we give some results 
that we will use in our analyses. For readers interested in the details of this formalism there are several good references 
including Refs. .1415, 16, 17. 

For a particle propagating from the initial point {xi,ti) to the final point (xf,tf) the transition amplitude is given 
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by the integral over all possible paths connecting the initial and the final points: 

where S and C denote, respectively, the classical action and the Lagrangian of the particle. The transition amplitude 
K{xf,tf;xi,ti) is called the propagator. It can be shown that for a general quadratic Lagrangian the form of the 
propagator reduces toi^ 

K{xf,tf-x,,U) = (2) 

where the factor (f){tf,ti) is a function of the initial and the final time. The subscript "cl" refers to classical solution 
of the equation of motion. 

An important property of the propagator, known as transitivity, is obtained by considering an instant of time t 
such that tf > t > ti 

K{xf,tf]Xi,ti)^ J dx K{xf,tf;x,t)K{x,t;Xi,ti). (3) 
We will make use of Eqs. ^ and ^ in our subsequent discussion. 



III. PATH INTEGRAL FORMULATION OF LINEARLY DAMPED SYSTEMS 



We now illustrate the method of path integrals outlined in Sec. II by applying it to some simple damped mechanical 
systems. We will assume that the motion takes place between fixed initial and final points and calculate the kernel 
for the systems. The kernel will then be used to study the evolution of a Gaussian wavepacket. 



A. Linearly Damped Harmonic Oscillator 



The equation of motion of a linearly damped harmonic oscillator is mx + l3x + muiQ^x — 0, where m is the mass of 
the particle, /3 is the damping coefficient, and loq is the frequency of its oscillations for P — 0. The general solution of 
the equation of motion for the over-damped (OD), critically damped (CD), and under-damped (UD) cases are shown 
in TableHl where A = /3/m, [xx better to write A instead xx] 7^ = — wq^) = — w^, and T = tf — ti. We use the 
boundary conditions x{ti) — x.^ and x{tf) = Xf to evaluate the integration constants A and B (see Table I). 

To construct the kernel we first need to know the Lagrangian, which is given by: 
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Note that the Lagrangian is explicitly time dependent. There are ways of choosing new coordinates so that the 
Lagrangian in Eq. ^ becomes time-independent It is easy to check that this Lagrangian reproduces the correct 
classical equation of motion for the damped harmonic oscillator. Several authors have looked at the path integral 
kernel for this Lagrangian. ^^'^^ A reasonably up-to-date review covering various aspects is available in Ref. [U We 
now write down the kernel for the under-damped case and then investigate the wavepacket evolution. The results for 
the other two cases are given in Tables |ll] and IIIII 

The classical action is evaluated by substituting the solution for the under-damped case given in Table [T] into Eq. ^ 
and integrating it over the time interval [ti,tf]. The result is 
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The effect of damping appears in Eq. ([5]) through the presence of A. In particular, the second term is entirely due to 
damping effects. It is interesting that if we use scaled coordinates Xi — x^e^*'/^ and Xf — XfC-^^f^^, we can rewrite 
the first term as a purely simple harmonic oscillator contribution. 

Because the Lagrangian is quadratic, the kernel is of the form given in Eq. ([2]). We make use of the transitivity of 
the kernel, that is, Eq. ([3|) to calculate (f>{tf,ti). After some algebra, we find: 



(l>{tf,t,) = ^itf,t)4>{t,t,) 
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which leads to 



The complete kernel turns out to be 
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where S is given by Eq. ([5]). 

To determine how a Gaussian wavepacket evolves for this kernel, we begin with the initial (t^ = 0) profile of the 
packet: 
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where tro^ is the variance of the Gaussian wavepacket, which is a measure of its width. Without any loss of generality 
we choose the wavepacket to be peaked at = a aX ti = 0. The wavepacket at a later time t is related to the 
wavepacket at ti = by 



K{xf, t; Xi,0)ip{xi, 0) dxi 



After some simplifications, we find 
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From Eq. (jlip we see that at any time t the wavepacket is peaked at 
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The wavepacket evolution is shown in Fig. [Tl and the dependence of the standard deviation at on t and A is shown 
in Fig. [2] From Fig. [Tl we notice that the width of the wavepacket pulsates and at various times it becomes less than 
the initial value ctq. From Fig. [5] we see that at shows the same behavior and after t ~ 11 remains less than ctq for 
A — 0.2 and wq = 0.5 (bold line). We see similar behavior for at for different values of A. The oscillations are less 
prominent for higher values of A. This behavior is seen for the critically and over-damped cases. From theoretical 
considerations, we expect that the under-damped case exhibits less prominent oscillations as A ^ 2u)q. We also note 
that for K — Q, at oscillates between ctq and (Tmax > (as is well known), but as the damping coefficient becomes 
nonzero, the variance drops below ctq a-t some time and goes to zero. This behavior coincides with the wavepacket's 
peak tending towards a; = 0. Thus, we conclude that the damping leads to localization of the particle around the 
minimum of the potential at a; = 0. 

The expectation value of x is 
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This result is the same as Eq. ([T^ and x{t) for the UD case in Table U if A and B are evaluated using the initial 
conditions x{Q) — a and i(0) — 0. That is, the peak of the wavepacket (corresponding to the maximum probability 
of finding the particle) follows the classical trajectory as expected. 



B. Uniform Gravitational Field with Linear Damping 



Consider a particle of mass m in a uniform gravitational field with a damping force proportional to its speed. This 
damping is an example of Stokes' law. The equation of motion of the particle is mx + px = mg, where (3 is the 
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damping coefficient and g the acceleration due to gravity. Recall that there is a terminal velocity, which the particle 
attains asymptotically. The general solution of the equation of motion is 



x{t) = A + Be 



-At 



^ A ' 
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where A — P/m^ and A and B arc integration constants. For the initial and final conditions x{ti) — Xi and x{tf) = xj, 
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The equation of motion (fTS]) can be derived from the Lagrangian 
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For this Lagrangian, the classical action in the time interval [U, tj] is 
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The calculation of the kernel can be done in a way similar to the the damped harmonic oscillator. We obtain 



™g r^-At. 



-At, 



K{tf,Xf;ti,Xi) 



TOAe^(*'+*/) 



2Trih{e 



Atf 



cAti 



iS/h 



(16) 



(17) 



(18) 



where S is given by Eq. p?)) . 

We will now consider the evolution of the wavepacket given in Eq. ^ . We make use of Eq. ITUl) and obtain 
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From Eq. (|19p we see that the wavepacket is peaked at 
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The wavepacket evolution is shown in Fig. [31 and the variation of at on the damping coefficient A and the time t 
is shown in Fig. U) From Fig. [3] we see that the width of the wavepacket increases initially and then becomes almost 
constant as the exponential part of Eq. ((20)) decays. From Fig. |4]we see that the variance at exhibits the same generic 
behavior for all values of A. However, the time taken to reach a near-constant value of a is different for different A 
values and represents the time needed to reach the terminal velocities in the corresponding cases. This behavior of 
the wavepacket seems to be characteristic of systems involving a terminal velocity though, as we show later there are 
interesting differences in the case for quadratic damping. The expectation value of x is 
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This result is the same as Eqs. (j2ip and (fT^l) if the constants are evaluated using the initial conditions a;(0) = a and 
i(0) = 0. 



IV. QUADRATIC DAMPING 



A. Path Integral Kernel and wavepacket Evolution 



We consider a particle moving in a uniform gravitational field with a damping force proportional to the square of its 
speed. Much work has been done on the quantization of this and similar systems with quadratic dampingj ^^i^^i^'^i^^i^^i^^ 
The equation of motion with this type of damping is mi-h/Ji^ = mg. Note that the equation is time-reversal invariant. 
The general solution is 



x{t) = [ln[cosh( V^Ai -I- A)]] B, 



(23) 
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where h. = (3/m and A and B are integration constants. The choice of a suitable Lagrangian for this case is interesting 
because there are nonequivalent Lagrangians which give rise to the same equation of motion. Consider for example 
the forms 
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To quantize the system we must judiciously choose the form that can be handled easily despite the fact that different 
Lagrangians can give rise to nonequivalent quantizations. The Lagrangian in Eq. (|25p is not so easy to use because of 
the presence of the square root in the path integral method. Thus, we choose the Lagrangian in Eq. (p4|) . Despite the 
presence of damping, the Lagrangians are not explicitly time-dependent unlike the damped harmonic oscillator or a 
particle in a gravitational field with linear damping. The Hamiltonian derived from Eq. I|24p is a conserved quantity, 
but does not correspond to the energy of the system. A discussion on the conserved quantities in damped systems is 
given in Ref. Is^ . 

We note that although the Lagrangian is not a quadratic Lagrangian, we can make it so by using the trans- 



formation: X = J e dx = e /A.^^ This transformation converts it into a Lagrangian similar to that of an simple 
harmonic oscillator with imaginary frequency whose results are knowni^ii^ or can be deduced from those of Sec. lIII 
by setting A = 0. 

We can write the action in terms of X and t as 



gAX^)dt. 



If we compare Eq. 



with the harmonic oscillator action given by: 
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we obtain uj = iyfgK — 17. We use the known results for the propagator of the harmonic oscillator and obtain the 
kernel in terms of X and 7 as 
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The evolution of a Gaussian wavepacket in the A-coordinate will be similar to Eq. 
A = 0). Therefore, in terms of the ^-coordinate we can write 
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If we set the exponent in Eq. 
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D) to zero, we see that the wavepacket is peaked at: 
Xf ^ ^ ^ ~^ ln(cosh7f). 



and Eq. (fTT|) (after setting 
(30) 

(31) 
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The wavepacket evolution is shown in Fig. [51 and the variation of Ut with time t and A is shown in Fig. [B) From 
Figs. [5] and [5] we see that the width of the wavepacket increases indefinitely and rapidly. From Fig. [S] we see that the 
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standard deviation at exhibits similar behavior for all values of the damping coefficient A. The only difference is the 
rate at which at grows, which can be derived from Eq. (j32p . Thus for motion under gravity, the linear and quadratic 
damping cases can be distinguished from each other by following the corresponding wavepacket evolution (see Fig. |4] 
for linear damping and Fig. [5] for quadratic damping). 

As before, we would like to calculate the expectation value of x. We first calculate (AT), which can be determined 
from Eq. (|14|) by setting A = 0. Using the relation X — exp^^ /A we get 

(e^^) = e'^"cosh7t. (34) 

If we expand the exponential on both sides and compare the coefficients of A, we obtain 

{x) ~ + ^ ln(cosh7i). (35) 

This result is the same as Eqs. ([55)) and (^5)). 



B. Connection with Field Theory 



In the following we will show that the classical equation of motion for a particle subject to damping proportional to 
the square of the velocity can be obtained from a zero-dimensional version of the field theory of tachyon matter that 
emerges from string theory^. This example provides a link between a field theory and a damped mechanical system. 
Recall other such connections such as that between the simple harmonic oscillator and the massive Klein-Gordon field 
theory. 

The action for the tachyon^^ matter field in a p + 1 dimensional spacetime is given as^ 

S = -j dP+\ V{T) ^l + rj'3d,TdjT, (36) 

where is the metric for a (p+l)-dimensional fiat spacetime with components rjij = diag(— 1, 1, . . .) for i, j — 0, . . . ,p; 
T(x) is the tachyon field, and V{T) ^ e"""^/'^ denotes the corresponding field potential. The value of the parameter 
a depends on the particular type of string theory of interest. 

We now consider the form of the action for a zero-space dimensional case, that is, p = 0. We identify the tachyon 
field with the coordinate x {T ^ x) in the classical problem. Note that in classical mechanics the action has the form 
S = jLdt and the zero-dimensional action has a similar form 

S^~jdt e-"^\/F^, (37) 

where the factor of 2 in the potential V{x) has been absorbed in a. The equation of motion that results from the 
action in Eq. (|57|) is 

X + ax^ = a. (38) 

We now scale x ^ bx and compare Eq. ([38)1 with the equation of motion x + Ax'^ = g. We obtain ab = A and a/b = g, 
which can be solved to yield a = ^fgA and b = ^ Aj g. We substitute these results into Eq. ([57]) and recover the 
Lagrangian in Eq. ([25)1 . Thus, we obtain a quadratically damped mechanical system out of a field theory. 



C. Damping as Geodesic Motion 



Another way of looking at the problem of quadratically damped motion is to picture it as the motion of a particle 
along a geodesic in a fictitious two-dimensional space. Consider the following form of a two-dimensional distance 
function (line element) 

ds2 ^ j(x)dx^ -t- h{x)dy^, (39) 

where y denotes the fictitious dimension and ]{x) and h(x) are unknown functions. Our aim is to show that for an 
appropriate choice of /(x) and /i(x), the equation of motion for a quadratically damped system can be identified as 
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a geodesic equation. For the metric 5^1/ = diag[f{x),h{x)], we calculate the nonzero components of the ChristofFel 
symbol^: 

f h' h' 

= — = r*' = = r'' f40'i 

XX yy 2/' 2h y^' ^ ' 

where the prime denotes the derivative of the function with respect to x. If we substitute these components into the 
well known geodesic equation given as 

(here A is any parameter on the geodesic, which in our case is the time t), we obtain the equation of motion along the 
two directions: 

y+^xy^O. (43) 
h 



If we integrate Eq. once, we have 



where C is an integration constant. We next substitute y in Eq. ([42]) and obtain 



which is the same equation as the quadratically damped equation of motion along the x direction provided that 
/'/2/ — A and C^h' /2fh? = g. These equations can be solved to reveal the form of the two functions: 

fix) = e^A- (46) 

and 

/i(a;) = -^e-2A-. (47) 
9 

It is now straightforward to calculate the components of the Riemann tensor, R°' p^vr^ Ricci tensor, i?^^ ~ R" nav, 
and the Ricci scalar, R = g>^'^R^^,-M 

9(^2 A3 

i?,, = -2A2, i?^^ ^ ^l_^^e-2A^ R = ~AA^e-^''\ (48) 

The presence of A in the curvature scalar R implies that the damping can be viewed as a curvature effect in this 
fictitious two-dimensional space. What distinguishes this case from the motion with linear damping discussed in 
section IIIB is that we cannot cast the equation of motion of the linear damping case in the form of a geodesic 
equation due to the absence of a term. In this sense, quadratic damping is unique. Thus the above connection 
provides us with additional geometric insight into the nature of quadratically damped motion. 



V. CONCLUDING REMARKS 



We have shown how to construct kernels for several damped mechanical systems and studied the evolution of a 
Gaussian wavepacket in each case. We demonstrated that for the linearly damped harmonic oscillator and a particle in 
a uniform gravitational field with linear and quadratic damping, we can see characteristic features of the damping from 
snapshots of wavepacket evolution. To motivate our consideration of quadratic damping, we related the corresponding 
equation of motion to a field theory and to geodesic motion in a fictitious two-dimensional space. 

Is a quadratically damped system damped? The Lagrangian is time independent and the system is Hamiltonian 
and conservative in the usual sense. The evolution of the wavepacket shows spreading, much like that of a free particle 
and unlike the linear damping system. We also note a similarity with the linearly damped system because in both 
cases, the particle attains a terminal speed. These issues suggest that it would be better to view the quadratically 
damped system as special and unlike the linearly damped case. 
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The components of the ChristofFcl symbol depend on the metric tensor components as F"^ = ^g"" [dp^Qav +dvgcrij. —dagfiv)- 
The components of Riemaim-Christoffel curvature tensor are given by R" = df^Tg^ — d^Tg^ + T'^^Tg^ — T^^Tg^. For 
details of its properties, see any text on the general theory of relativity. 
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Tables 



Case 


x(t) 




B 


OD 


g- At/2 gj^g].^ ^ g 


[x^ e^*i sinh ft j —x j- e"^*/ ^ sinh 7t j] 
sinh^T 


[x/e^*//^ cosh -fti - Xie^'^l^ cosh ^t,] 


CD 


(-''-'/''[A + Bt] 






UD 


g-Ai/2 u)t + B sin tjt] 


r At,-/2 - i A.tf/-2 . , 
[XjC *' sinwty— xye J sinujtj] 

sinuT 


sinl.T [^/fi"^'^^^ COStjti - l-je'^**''^ COSt^t/] 



TABLE I: The general solutions of the equations of motion of a damped harmonic oscillator in different cases. We show the 
expressions for x{t), A, and B for the different cases. (See the discussion in section IIIA.) 



Case 


S (Action) 


K (Kernel) 


CD 


^[xh^'^ + x}e^'f - 2x,Xfe'"^'^+'f ^/^] + ^(xh^'^ - x}e^'s) 


V 2TTihT 


OD 


^^[(4e^*- +x^e'^*/)cosh7T- 2x,x,e^(*.+V)/2]+ 


1 A(tj+tf)/2 i „ 
V 27riBsinh7T 



TABLE II: Action (5*) and Kernel {K) for the critically and over-damped harmonic oscillator. 



Case 



CD 



:exp[ 



Xf—ae 



OD 



2a 



ao^e-'^*[(cosh7t + ^sinh7t)^ + (li^)'] 



TABLE III: Probability amplitude \il}{xf,t)f and measure of dispersion for the critically (CD) and over-damped (OD) 
harmonic oscillator. 
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Figure Captions 



Wavepacket Evolution for Under-Damped Harmonic Oscillator 
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FIG. 1: Evolution of wavepacket in a liarmonic oscillator potential with linear damping for the parameters A = 0.2, uo = 0.5, 
and a = 1.0. 
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Wavepacket Evolution for Linearly Damped Motion under Gravity 



70 - 



60 ■ 



50 



t=0 
0=0 



t 

40 ■ 



30 ■ 



)050 



20 - 



10 ■ 



t=0.4 
0=0.0063 




-0.05 0.05 0.6891 



1.926 4.242 
X -> 



t=7 

0=0.0126 




t=13 

0=0.0126 




t=20 

0=0.0126 




58.81 



117.6 186.2 



FIG. 3: Evolution of wavepacket for motion under gravity with linear damping for the parameters A — 1.0, a — 0.0, and 
5 = 9.8. 
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Evolution of Std. Deviation w.r.t.A for Linearly Damped Motion under Gravity 




2 




FIG. 4: Variation of at with t and A for motion under gravity with linear damping [Eq. (|20|) ]. The value of A for each curve is 
shown just above the corresponding curve. 



80 r 



70 ■ 



60 



50 - 



T 

cv!_ 40 ■ 



30 



20 - 



10 ■ 



OL 



t=0 
a=0 



1)050 



Wavepacket Evolution for Quadratically Damped Motion under Gravity 
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FIG. 5: wavepacket evolution for motion under gravity with quadratic damping with A — 1.0, a = 0.0, and g = 9.1 




FIG. 6: Variation of at with t and A for motion under gravity with quadratic damping, [Eq. (|32[)]. The value of A for each 
curve is shown just above the corresponding curve. 



